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I.  Research  Activities 


As  listed  below,  we  have  completed  all  our  goals  in  our  research  toward  the  "Multiscale, 

Intermittent,  Turbulent  Fluctuations  in  Space  Plasmas  and  Their  Influence  on  the  Interscale 

Behavior  of  the  Space  Environment 

•  We  developed  an  innovative  method  of  analysis  for  intennittent  fluctuations  called  ROMA 
(rank-ordered  multifractal  analysis).  The  method  allows  us  to  calculate  the  implicit 
multifractal  spectrum  and  crossover  functions  for  fully  developed  intennittent  turbulence  in 
space  plasmas. 

•  We  applied  the  ROMA  method  to  situations  where  the  physical  domain  changes  from  the 
kinetic  size,  to  intermediate  ranges,  and  all  the  way  to  the  macroscopic  magnetohydrodynamic 
region.  We  developed  global  scale  invariant  theories  to  map  the  probability  distribution 
functions  of  fluctuations  of  all  sizes  onto  one  master  scaled  distribution. 

•  We  performed  detailed  analyses  of  magnetic  field  intensity  fluctuations  in  the  magnetospheric 
cusp  and  sheath  regions  based  on  data  obtained  by  the  Cluster  instruments.  The  results 
revealed  interesting  inhomogeneous  intermittent  turbulence  in  nonuniform  magnetized 
environments  as  predicted  by  ROMA. 

•  We  studied  the  intennittent  turbulence  in  the  auroral  zone  based  on  actual  rocket  data  kindly 
provided  to  us  by  our  colleagues  (Paul  Kintner  and  Eric  Klatt)  at  the  Cornell  University.  We 
also  studied  and  compared  the  multifractal  nature  of  both  the  results  from  our  direct  numerical 
simulations  and  in-situ  rocket  observations  using  ROMA. 

•  We  developed  a  new  imaginative  theory  called  DSRG  (direct  simulation  renormalization 
group).  It  incorporates  large  scale  numerical  simulation  with  the  most  sophisticated  theory  in 
statistical  physics,  the  dynamic  renormalization  group.  Using  DSRG,  we  were  able  to 
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calculate  and  predict  complicated  multifractal  properties  of  intennittent  turbulence  properties 
of  space  plasmas  under  various  realistic  conditions. 

•  We  made  definitive  strides  toward  the  basic  understanding  of  coarse-grained  dissipation 
process  that  can  lead  to  the  merging  of  coherent  plasma  structures  resulting  in  the 
development  of  intermittent  turbulence  as  well  as  large-scale  magnetic  reconfigurations  such 
as  reconnections  in  the  Earth's  magnetotail. 

•  We  obtained  large-scale  numerical  simulation  results  of  the  fluctuation  probability 
distribution  functions  and,  in  particular,  the  local  intennittency  measures  based  on  the  wavelet 
transfonns  for  intermittent  turbulence.  These  results  are  obtained  using  the  concept  of 
interacting  plasma  coherent  structures.  The  simulation  results  demonstrate  clearly  the 
sporadic,  localized,  and  scale-dependent  nature  of  intermittent  turbulence  in  space  plasmas. 

•  The  concepts  of  scale  invariance  and  symmetry-breaking  relating  to  the  basic  properties  of 
intennittent  turbulence  are  studied  based  on  the  dynamical  theory  of  the  renormalization 
group.  Our  results  from  ROMA  and  DSRG  convincingly  demonstrated  that  the  phenomena 
of  self-organized  criticality  and  intermittent  turbulence  are  intimately  related  through  the 
concept  of  crossover  processes  of  scale  invariants. 

Throughout  the  period,  the  MIT  group  actively  interacted  with  Dr.  J.  R.  Jasperse  and  his 
theoretical  space  plasma  group  at  the  Air  Force  Research  Laboratory  at  Hanscom  Field  in 
Massachusetts  (now  at  Boston  College).  This  joint  collaborative  effort  led  to  the  important 
understanding  of  the  process  of  particle  interactions  with  intennittent  plasma  turbulence, 
particularly  in  the  return  current  region  of  the  auroral  zone.  This  beneficial  collaboration 
continued  throughout  the  research  grant  period. 
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Such  efforts  as  listed  above  have  resulted  in  a  number  of  peer-reviewed  research 
publications  as  detailed  in  Section  II.  In  addition,  we  have  been  invited  to  deliver  a  number  of 
specialty  lectures  on  these  subjects  at  major  institutions  and  international  conferences. 

A  "brief  tutorial  review"  recently  published  in  Nonlinear  Processes  in  Geophysics, 
"Rank-ordered  Multifractal  Spectrum  for  Intennittent  Fluctuations",  is  enclosed  to  highlight  our 
major  activities  and  findings  sponsored  by  this  research  grant. 
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Abstract.  Intermittent  fluctuations  are  the  consequence  of 
the  dynamic  interactions  of  multiple  coherent  or  pseudo- 
coherent  structures  of  varied  sizes  in  the  stochastic  me¬ 
dia  (Chang,  1999).  We  briefly  review  here  a  recently  de¬ 
veloped  technique,  the  Rank-Ordered  Multifractal  Analysis 
(ROMA),  which  is  both  physically  explicable  and  quantita¬ 
tively  accurate  in  deciphering  the  multifractal  characteristics 
of  such  intermittent  structures  (Chang  and  Wu,  2008). 

The  utility  of  the  method  is  demonstrated  using  results  ob¬ 
tained  from  large-scale  2-D  MHD  simulations  as  well  as  in- 
situ  observations  of  magnetic  field  fluctuations  from  the  in¬ 
terplanetary  and  magnetospheric  cusp  regions,  and  the  broad¬ 
band  electric  field  oscillations  from  the  auroral  zone. 


1  Introduction 

As  plasmas  interact  nonlinearly  in  the  electromagnetic  envi¬ 
ronment  of  space,  coherent  structures  are  generally  created. 
For  example,  in  the  solar  wind,  field-aligned  magnetic  flux 
tubes  of  varied  sizes  have  been  observed  (Bruno  et  al.,  2001). 
Other  examples  of  coherent  structures  observed  in  space  are 
double  layers,  ion  holes,  electron  holes,  drift  Alfven  vor¬ 
tices,  current  filaments,  stationary  and  convecting  structures. 
When  these  (sometimes  partially  formed)  structures  interact 
stochastically,  intermittent  fluctuations  are  generated. 
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For  intermittent  turbulence,  one  may  visualize  the  fluctu¬ 
ations  to  be  composed  of  many  types,  each  being  character¬ 
ized  by  a  particular  fractal  dimension.  Two  questions  arise: 

(i)  What  are  the  different  types  of  fractal  dimensions?  and 

(ii)  How  are  they  distributed  in  the  turbulent  medium?  Re¬ 
cently,  a  new  technique  of  analyzing  intermittent  fluctuations 
has  been  developed  to  specifically  address  these  questions 
(Chang  and  Wu,  2008).  The  technique,  known  as  Rank- 
Ordered  Multifractal  Analysis  (ROMA),  retains  the  spirit  of 
the  traditional  structure  function  analysis  and  combines  it 
with  the  idea  of  one-parameter  scaling  of  monofractals. 

We  briefly  review  this  method  below  in  terms  of  large- 
scale  two-dimensional  (2-D)  magnetohydrodynamic  (MHD) 
simulations  and  then  apply  it  to  the  study  of  several  observa¬ 
tional  examples  in  the  space  environment. 

2  Fractals 

We  begin  with  an  elemental  introduction  of  the  concept  of 
fractals.  In  ordinary  Euclidean  geometry,  the  area  of  a 
two-dimensional  square  may  be  expressed  as  the  square  of 
the  length  of  one  of  its  sides,  and  the  volume  of  a  three- 
dimensional  cube  is  the  cube  of  the  length  of  one  of  its  sides. 
Generalizing,  the  volume  V  of  a  hyper-dimensional  cube  of 
dimension  d  is  V  —  Ld ,  where  L  is  the  length  of  one  of 
the  sides  of  the  hypercube.  Consider  now  the  fluctuations 
S X  of  some  physical  quantity  A  of  a  turbulent  plasma.  We 
may  construct  some  convenient  scale-dependent  measure  p 
of  the  fluctuations  -  equivalent  to  the  volume  of  the  hyper¬ 
cube  -  and  study  how  the  measure  varies  with  the  scale  size 
S  -  equivalent  to  the  length  of  one  side  of  the  hypercube. 
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Fig.  1.  2-D  MHD  simulation  results  with  (512x512)  grid  points  for  the  PDFs  of  the  magnetic  field  energy  density  fluctuations  SB2  at  scales 
of  S  =  2  (red),  4  (green),  8  (black),  16  (blue)  and  32  (magenta)  units  of  grid  spacing.  Left  panel  displays  the  calculated  results  of  P(SB2,S ) 
and  the  right  panel  displays  the  corresponding  scaled  distributions  for  s  =  0.34. 


If  the  variation  between  /r  and  <5  follows  that  of  a  power  law, 
H  ~  8d{,  then  the  value  d{  may  be  visualized  as  the  dimension 
of  the  fluctuations  for  the  particular  chosen  measure  fi. 

If  the  obtained  dimension  of  the  measure  of  the  fluctua¬ 
tions  may  be  understood  in  terms  of  straight  forward  dimen¬ 
sional  analysis,  one  visualizes  the  fluctuations  fully  occupy¬ 
ing  the  space  in  which  they  reside.  On  the  other  hand,  if 
the  obtained  dimension  of  the  fluctuations  differs  from  that 
surmised  from  dimensional  analysis,  the  fluctuations  do  not 
fully  occupy  the  space  in  which  they  exist.  Such  a  dimension 
has  been  termed  a  “fractal”  or  fractal  dimension  by  Mandel¬ 
brot  (1977). 

We  generalize  this  idea  to  the  study  of  the  probability  dis¬ 
tribution  functions  (PDFs)  of  fluctuations.  To  be  more  spe¬ 
cific,  we  refer  to  the  results  of  a  2-D  MHD  simulation  for 
homogeneous  turbulence  (Chang  et  al.,  2004).  In  the  simu¬ 
lation,  ideal  compressible  MHD  equations  expressed  in  con¬ 
servative  forms  were  solved  numerically  in  a  doubly  periodic 
(x,z)  domain  of  length  2tt  in  both  directions  with  512  by 
512  grid  points  (Fig.  1)  or  1024  by  1024  grid  points  (Fig.  2) 
using  the  WENO  code  (Jiang  and  Wu,  1999)  so  that  the  total 
mass,  energy,  magnetic  fluxes  and  momenta  were  conserved. 
The  initial  conditions  consisted  of  random  magnetic  and  ve¬ 
locity  fields  with  a  constant  total  pressure  for  a  high  beta 
plasma.  (We  note  that  the  ideal  MHD  equations  are  scale 
independent.  They  can  be  made  dimensionless  with  the  spec¬ 
ification  of  the  units  of  length,  mass  and  time.  For  reference, 
the  initial  parameters  in  the  simulation  for  the  maximum  per¬ 
turbation  velocity  v  and  magnetic  field  B  were  chosen  such 
that  the  sound  wave  and  Alfven  wave  traveling  times  through 
a  distance  2 jt  are  about  4.4  and  60,  respectively.) 

After  sufficient  elapsed  time,  the  system  evolved  into  a  set 
of  randomly  interacting  multiscale  coherent  structures  ex¬ 
hibiting  classical  aspects  of  intermittent  fluctuations.  Spa¬ 
tial  values  of  the  magnetic  field  energy  density  Br  were 
collected  over  the  entire  (x,  z)-plane  at  a  given  time  and 
PDFs  P(8B2 ,8 )  were  constructed  for  the  spatial  fluctuations 
8B2  at  scale  intervals  5  =  2,4,8,16,32  grid  points.  Fig.  1 
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Fig.  2.  Rank-ordered  multifractal  spectrum  for  a  (1024  x  1024)  2-D 
MHD  simulation,  s(T),  Y  =  Ii5B2|/<5s;  |(5B2|  in  units  of  bin  size 
and  8  in  grid  spacing.  The  spectrum  is  calculated  for  8  contiguous 
ranges  of  AL. 

(left  panel).  The  PDFs  are  seen  to  be  non-Gaussian  and  be¬ 
come  more  and  more  heavy-tailed  at  smaller  and  smaller 
scales. 

We  now  attempt  to  understand  the  scaling  behavior  among 
the  family  of  PDFs  of  the  turbulent  fluctuations  in  analogy 
to  the  idea  of  fractals.  We  assume  that  the  8B2’s  and  the 
PDFs  vary  with  8  as  power  laws:  8B2/8a  —  I  and  P/8b  —  J, 
where  ( a,b )  are  the  exponents  and  (/,  J )  are  constants  -  i.e. 
invariants  with  respect  to  the  scale  8.  If,  in  addition,  the  form 
of  P(8B2,8 )  is  invariant  with  respect  to  the  changes  of  the 
scale  8,  then  there  must  be  a  functional  relation  between  the 
two  invariants,  (/,  J),  (Chang  et  al.,  1973).  Therefore,  along 
with  the  normalization  condition  for  the  PDFs,  we  obtain  the 
following  scaling  form  for  the  PDFs: 

(1) 

where  s  —  a  —  —b  is  the  lone  scaling  exponent,  and  Ps 
is  a  scaling  function  of  the  invariant  SB2  —  8B2 /8s  =  Y . 
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The  scaling  exponent  5  may  be  interpreted  as  the  fractal 
(monofractal)  measure  for  (1).  If  the  PDFs  are  self-similar 
Gaussian  distributions  for  all  scales  representing  random  dif¬ 
fusion,  the  scaling  exponent  s  is  known  to  be  equal  to  0.5. 
For  other  monofractal  distributions,  the  scaling  exponent 
may  take  on  any  positive  real  value. 

An  attempt  to  collapse  the  unsealed  PDFs  of  Fig.  1  (left 
panel)  according  to  expression  (1)  indicated  approximate 
scaling  with  an  estimated  scaling  exponent  s  —  0.34,  Fig.  1 
(right  panel).  The  scaling  was  quite  robust  for  small  values 
of  <5B2,  but  became  less  and  less  satisfactory  as  the  value  of 
<5B2  increased. 

When  the  scaling  law  (1)  is  fully  satisfied  for  the  full 
range  of  the  scale  invariant,  the  scaling  of  the  PDFs  is  called 
self-similar  or  monofractal.  Otherwise,  the  fluctuating  phe¬ 
nomenon  is  multifractal. 

3  Multifractals  and  ROMA 

One  established  procedure  for  measuring  the  multifractal  na¬ 
ture  of  fluctuations  is  by  studying  the  moments  of  the  PDFs, 
Sm  -  commonly  called  structure  functions.  If  the  scaling  ex¬ 
ponents  Km  for  Sm  vary  linearly  with  the  moment  order  m, 
one  visualizes  that  the  various  portions  of  the  PDFs  repre¬ 
senting  the  different  sizes  of  the  fluctuations  have  the  same 
fractal  characteristics  and  therefore  all  the  PDFs  collapse 
onto  one  scaled  PDF  according  to  (1).  In  fact,  one  may  easily 
demonstrate  that  PDFs  satisfying  (1)  will  have  the  monofrac¬ 
tal  property  Km  —  sm,  and  s  has  been  called  the  Hurst  expo¬ 
nent. 

For  multifractal  fluctuations,  the  fluctuations  within  a 
small  differential  range  dY  about  the  value  Y  may  be  as¬ 
sumed  to  exhibit  approximate  monofractal  behavior.  Thus, 
the  differential  structure  function  dSm  will  vary  with  the 
scale  as  Ssm  within  the  differential  range  according  to  (1), 
i.e., 

dSm  =  (sB2^"  p(sB2,S}d(sB2^)  =8smYm  Ps(Y)dY  (2) 

The  solution  s(Y)  for  the  above  functional  differential  ex¬ 
pression  (if  the  ansatz  is  valid)  may  be  obtained  approxi¬ 
mately  by  integrating  (2)  over  small  contiguous  ranges  of 
AT  with  the  assumption  that  within  each  incremental  range 
the  scaling  exponent  s  is  approximately  constant.  Figure  2 
shows  the  results  of  such  a  rank-ordered  spectmm  using  the 
numerical  data  from  a  MHD  simulation  similar  to  that  dis¬ 
cussed  in  the  previous  section.  Because  of  the  approximate 
symmetry  of  the  PDFs  with  respect  to  SB2,  we  have  defined 
Y  =  |<5B2|/<5T  We  note  from  Fig.  2  that  the  scaling  expo¬ 
nent  s  appears  to  approach  the  self-similar  value  of  0.5  rep¬ 
resenting  Gaussian  random  diffusion  as  Y  — »•  0.  As  the  size 
of  the  scaled  fluctuation  Y  increases  the  scaling  exponent  de¬ 
creases  indicating  that  the  fluctuations  are  becoming  less  and 
less  space-filling. 
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This  method  of  analyzing  intermittent  fluctuations, 
ROMA  (rank-ordered  multifractal  analysis),  was  first  sug¬ 
gested  by  Chang  and  Wu  (2008). 

4  Application  to  space  plasma  observations 

Recently,  ROMA  has  been  applied  to  the  analyses  of  turbu¬ 
lence  data  observed  within  several  regimes  of  the  space  en¬ 
vironment.  We  briefly  review  these  results  below. 

4.1  Solar  wind 

Most  of  the  turbulence  data  in  space  came  from  measure¬ 
ments  obtained  by  one  single  spacecraft.  For  solar  wind, 
the  collected  data  are  generally  analyzed  based  on  the  Tay¬ 
lor  hypothesis  (1938)  which  assumes  that  the  transit  time  of 
eddies  (bundled  fluctuations)  is  much  less  than  the  characte¬ 
ristic  evolution  time.  For  fully  developed  MHD  turbulence, 
the  characteristic  eddy  turnover  (evolution)  time  may  be  es¬ 
timated  as 

(X/2jt)/[(SB/B0)va] 

where  X  is  the  typical  size  of  the  eddy,  va  is  the  Alfven  speed 
based  on  the  mean  magnetic  field  Bo  and  SB  is  the  average 
fluemating  magnetic  field  of  the  eddy  in  question.  The  tran¬ 
sit  time  of  the  eddy  is  X/vs  where  us  is  the  relative  stream¬ 
ing  speed  between  the  plasma  medium  and  the  spacecraft. 
Thus,  as  it  has  been  demonstrated  by  Matthaeus  and  Gold¬ 
stein  (1982)  that  the  requirement  of  f (evolution)  be  much 
larger  than  /(transit)  becomes  v$/va  »  2tt  ( SB/Bq ).  This 
condition  is  generally  satisfied  for  solar  wind  turbulence. 

Because  the  spacecraft  speed  is  much  smaller  than  the  so¬ 
lar  wind  speed,  the  data  gathered  by  instruments  aboard  the 
spacecraft  essentially  sample  the  spatial  fluctuations  in  the 
direction  of  the  solar  wind,  particularly  at  times  when  the  so¬ 
lar  wind  turbulence  is  homogeneous  at  large  scales  and  fully 
developed.  One  may  then  construct,  for  example,  the  proba¬ 
bility  distribution  function  P(SB2,  r)  with  r  being  the  time 
scale  of  the  fluctuations  in  analogy  to  the  analysis  of  the  si¬ 
mulation  results  discussed  above. 

Based  on  a  particular  set  of  WIND  solar  wind  magnetic 
field  data,  Hnat  et  al.  (2002)  demonstated  that  the  PDFs  of  the 
magnetic  field  energy  density  fluctuations  P(SB2,t)  nearly 
collapsed  onto  a  single  scaling  curve  according  to  (1),  in¬ 
dicating  that  the  observed  solar  wind  turbulence  was  essen¬ 
tially  monofractal. 

Using  a  different  set  of  WIND  data,  Chang  et  al.  (2008) 
obtained  a  ROMA  spectrum  for  the  observed  solar  wind  tur¬ 
bulence,  Fig.  3.  The  spectrum  had  the  signature  of  that  of 
Fig.  2.  However,  it  was  much  flatter,  i.e.,  the  range  of  va¬ 
lues  of  the  fractal  index  s  (roughly  between  0.44  and  0.38) 
for  the  full  observed  range  of  the  scale  invariant  Y  =  SB2 /ts 
was  quite  small.  In  fact,  it  was  found  that  the  observed  PDFs 
nearly  collapse  unto  one  single  curve  according  to  (1)  by  set¬ 
ting  the  fractal  index  s  as  a  constant  (0.44),  Fig.  4,  quite  si¬ 
milar  to  the  result  obtained  by  Hnat  et  al.  (2002). 
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,  r  Fig.  4.  Scaled  PDFs  for  s  =  0.44  for  the  solar  wind  data  of  Fig.  3. 

Fig.  3.  Rank-ordered  multifractal  spectrum  for  a  particular  set  of 

solar  wind  data  from  WIND  for  fluctuations  between  3  and  1500  s. 

The  spectrum  is  calculated  for  5  contiguous  ranges  of  AF. 


We  now  proceed  to  consider  turbulent  fluctuations  within 
other  space  environments  where  multifractal  structures  are 
more  evident. 

4.2  Magnetospheric  cusp  turbulence 

One  of  its  regions  of  the  magnetosphere  most  exposed  to  the 
solar  wind  is  the  magnetospheric  cusp.  It  has  been  demon¬ 
strated  by  Echim  et  al.  (2007)  that  in  the  cusp  the  magnetic 
energy  density  fluctuations  have  non-Gaussian  PDFs  with 
large  extended  tails  at  all  scales  (time  lags).  Fig.  5.  And  the 
PDFs  do  not  scale  with  respect  to  their  variances,  a  strong 
indicator  of  intermittency  -  i.e.,  not  self-similar  under  scale 
transformations. 

There  was  strong  indication  that  the  statistical  properties 
of  the  fluctuations  in  the  cusp  are  inhomogeneous  and  did  not 
represent  fully  developed  turbulence.  How  inhomogeneity 
may  be  included  explicitly  and  more  rigorously  in  a  ROMA 
analysis  is  discussed  in  the  concluding  section. 

A  demonstration  of  the  validity  of  the  Taylor  hypothe¬ 
sis  is  not  straightforward  for  the  cusp  region.  In  Echim  et 
al.  (2007),  we  provided  an  argument  for  the  justification  of 
the  hypothesis  to  the  data  set  that  was  analyzed.  Firstly, 
Lundin  et  al.  (2003)  found  that  the  plasma  bulk  velocity  is 
of  the  order  of  1 00-200  km/s  for  the  same  Cluster  data  as 
those  discussed  in  this  paper.  Thus  the  plasma  speed  took  on 
values  larger  than  the  Cluster  velocity  indicating  that  tem¬ 
poral  variations  could  be  transformed  into  spatial  variations. 
Secondly,  a  correlation  analysis  was  performed.  Both  the 
dynamic  linear  autocorrelation  factors  for  the  B-field  fluctu¬ 
ations  on-board  each  Cluster  satellite  and  the  dynamic  cross 


correlation  factors  for  the  same  B-field  component  (or  am¬ 
plitude  or  B2)  between  pairs  of  Cluster  spacecraft  were  com¬ 
puted.  It  was  found  that  the  dynamic  autocorrelation  coeffi¬ 
cients  are  similar  to  the  cross  correlation  factors,  which  could 
be  interpreted  as  a  verification  that  temporal  fluctuations  cor¬ 
respond  to  spatial  fluctuations  (Weygand  et  al.,  2005). 

With  this  assumption,  a  ROMA  analyses  had  been  per¬ 
formed  for  the  chosen  data  set  in  the  cusp  region.  Fig.  6 
(Famy  et  al.,  2008).  When  the  turbulence  is  fully  devel¬ 
oped,  the  ROMA  spectrum  generally  has  a  smooth  mono¬ 
tonic  behavior.  On  the  other  hand,  if  the  turbulence  is  un¬ 
stable  and  not  yet  fully  developed,  its  ROMA  spectrum  may 
vary  rapidly  over  small  ranges  of  T.  Under  such  circum¬ 
stances,  the  increment  AT  in  the  ROMA  analysis  need  to  be 
adjusted  sufficiently  small  to  avoid  the  spurious  occurrence 
of  double  roots.  Since  such  circumstances  occur  generally 
for  small  Y  values,  there  usually  will  be  sufficient  number 
of  data  points  to  determine  the  spectrum  values  even  for  rel¬ 
atively  small  AT  increments  (Tam  et  al.,  2010).  We  have 
employed  this  procedure  of  variable  AT  in  determining  the 
ROMA  spectra  in  this  and  the  following  section. 

We  note  that  for  small  values  of  T,  the  fluctuations  were 
persistent  ( s  >0.5),  indicating  the  turbulence  was  unstable 
and  probably  not  yet  completely  fully  developed.  For  larger 
values  of  T,  the  fluctuations  became  anti-persistent  ( s  <  0.5) 
and  the  turbulence  was  probably  well-developed  and  be¬ 
came  sparser  and  sparser  as  the  value  of  T  increased  simi¬ 
lar  to  the  MHD  simulation  result.  The  ROMA  spectra  for 
all  4  spacecraft  were  very  similar;  indicating  the  turbulence 
there  was  essentially  isotropic. 
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<5b2=(AB2-<AB2>)/(i  (5b2=(AB2— <AB2>)  /  a 


Fig.  5.  Typical  PDFs  of  the  magnetic  energy  density  fluctuations  measured  by  the  four  Cluster  satellites  in  the  cusp;  the  PDFs  have  been 
scaled  with  respect  to  their  variance.  Dashed  curves  are  Gaussian  distributions.  The  scales  (r)  are  color  coded  and  given  in  seconds. 


4.3  Auroral  electric  field  fluctuations 

Tam  et  al.  (2010)  have  applied  the  ROMA  technique  to 
analyze  the  multifractal  characteristics  of  the  auroral  zone 
electric  field  fluctuations  observed  by  the  SIERRA  (Sound¬ 
ing  of  the  Ion  Energization  Region:  Resolving  Ambiguities) 
rocket.  Since  the  dominant  fraction  of  the  electric  field  vari¬ 
ations  in  the  observational  region  were  known  to  be  nearly 
non-propagating  intermittent  turbulence  composed  mainly  of 
multiple  interacting  inertial  Alfvenic  potential  structures,  the 
observed  time  series  were  interpreted  in  terms  of  spatial  fluc¬ 
tuations  with  A«1/t,  where  f/  ~  1 .5  km/s  is  the  horizontal 
speed  of  the  rocket,  r  the  time  scale  and  A  the  corresponding 
spatial  scale. 

In  analyzing  the  data,  it  was  discovered  that  the  observed 
fluctuations  (commonly  called  the  broad  band  extremely  low 
frequency  or  BBELF  fluctuations)  span  across  four  contigu¬ 
ous  regimes  of  scales  (from  kinetic  to  inertial  Alfvenic,  to 
intermediate,  and  eventually  to  the  MHD  scales)  with  dif¬ 
ferent  multifractal  characteristics.  A  ROMA  spectrum  was 
found  for  each  regime,  5,  (7,)  for  i  —  1,2, 3, 4.  The  dividing 
times  between  the  adjacent  regimes  were  determined  to  be 
roughly  r  ~  80  ms,  r  ~  160  ms,  and  r  bs  320  ms. 


CLUSTER  1 
CLUSTER  2 
CLUSTER  3 
CLUSTER  4 
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Fig.  6.  Rank-ordered  multifractal  spectra  for  a  typical  Cluster  pass 
in  the  cusp. 
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Fig.  7.  Top:  rank-ordered  spectra  .sq  (black),  S2  (red),  S3  (blue), 
and  S4  (green)  over  the  ranges  of  Tgiobal-  Bottom:  global  scaling 
function  Ps j  (Tgiobal)  obtained  from  the  PDFs  at  all  the  time  scales 
of  the  four  regimes. 

Thus,  the  observed  intermittent  turbulence  of  the  electric 
field  fluctuations  were  found  to  be  characterized  by  two  rank- 
ordered  parameters:  the  index  i  and  the  power-law  scaling 
variable  F,- .  Assuming  the  nonlinear  crossover  regions  be¬ 
tween  the  adjacent  regimes  were  narrow,  the  ROMA  spectra 
were  mapped  onto  one  global  scaling  variable  with  a  single 
global  scaling  function.  Fig.  7.  Such  a  mapping  served  the 
purpose  of  demonstrating  the  ROMA  spectrum  as  multival¬ 
ued  and  provided  a  convenient  display  of  the  individual  mul¬ 
tifractal  characteristics  for  each  rank-ordered  index  i.  It  also 
demonstrated  the  generalization  of  the  idea  of  monofractal 
mapping  (1)  to  the  concept  of  multi -parameter  mapping  in¬ 
volving  crossover  effects  among  multifractals. 

5  Summary,  conclusion  and  speculative  remarks 

1 .  We  briefly  reviewed  the  ROMA  technique  of  analyzing 
multifractals.  The  method  was  developed  using  the  re¬ 
sults  of  large  scale  2-D  MHD  simulations.  Applications 
of  the  method  to  the  interplanetary  medium,  the  magne- 
tospheric  cusp,  and  the  auroral  region  were  briefly  de¬ 
scribed. 
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2.  Such  an  implicit  multifractal  spectrum  has  several  ad¬ 
vantages  over  the  results  obtainable  using  the  conven¬ 
tional  structure  functions.  Firstly,  the  utility  of  the  spec¬ 
trum  is  to  fully  collapse  the  unsealed  PDFs.  Secondly, 
the  physical  interpretation  is  clear.  It  indicates  how  in¬ 
termittent  (in  terms  of  the  value  of  s )  the  scaled  fluc¬ 
tuations  are  once  the  value  of  F  is  given.  Thirdly, 
the  determination  of  the  values  of  the  fractal  nature  of 
the  grouped  fluctuations  is  not  affected  by  the  statistics 
of  other  fluctuations  that  do  not  exhibit  the  same  frac¬ 
tal  characteristics.  Fourthly,  it  provides  a  natural  con¬ 
nection  between  the  one-parameter  monofractal  scaling 
idea  (1)  and  the  multifractal  phenomenon  of  intermit¬ 
tency. 

3.  Most  of  the  existing  analyses  on  multifractals  are  in 
terms  of  either  spatial  or  temporal  fluctuations.  On 
the  other  hand,  fluctuations  that  exist  in  the  real  world 
are  generally  spatiotemporal  in  nature.  In  an  interest¬ 
ing  paper,  Pulkkinen  et  al.  (2006)  discussed  the  scaling 
of  the  second  order  structure  functions  that  are  intrin¬ 
sically  spatiotemporal  and  monofractal  or  self-affine. 
In  principle,  ROMA  may  be  generalized  to  consider 
the  scaling  of  multifractals  in  a  similar  manner.  For 
monofractal  processes,  the  underlying  stochastic  equa¬ 
tion  may  generally  be  expressed  in  terms  of  a  Langevin 
equation  with  Gaussian  noise,  or  equivalently  a  Fokker- 
Planck  equation.  For  multifractal  fluctuations,  the  un¬ 
derlying  stochastic  equation  would  be  a  spatiotemporal 
Langevin-like  equation  driven  by  a  generalized  noise. 
And  the  solution  of  such  an  equation  would  require 
the  meticulous  application  of  the  methods  of  path  inte¬ 
grals  and  the  nonlinear  dynamic  renormalization-group 
(Chang  et  al.,  1978,  1992,  2004).  The  results  of  a  multi¬ 
parameter  ROMA  study  with  nonlinear  crossover  ef¬ 
fects  based  on  the  concept  of  spatiotemporal  fluctua¬ 
tions  can  shed  light  to  the  description  of  such  a  theoret¬ 
ical  analysis.  These  studies  are  currently  being  pursued 
by  the  authors. 

4.  Traditional  multifractal  analyses  are  based  on  the  statis¬ 
tics  of  incremental  changes.  An  alternative  popular 
method  of  complexity  analysis  of  fluctuations,  the  con¬ 
cept  of  self-organized  criticality  (SOC)  (Bak  et  al., 
1987),  considers  the  statistical  description  of  avalanches 
structured  around  the  ideas  of  threshold  cutoffs  and 
finite-size  scaling  effects.  Recently,  Uritsky  et  al.  (2007, 
2009a)  analyzed  an  extended  data  set  of  extreme  ultra¬ 
violet  images  of  the  solar  corona  using  both  methods. 
The  data  exhibited  SOC  power-law  avalanche  statistics 
as  well  as  multifractal  properties  of  structure  functions 
for  spatial  activity.  It  will  be  interesting  to  investigate 
if  this  simultaneous  existence  of  SOC  and  multifrac¬ 
tal  phenomenon  is  universal  for  dynamical  complexity 
processes  or  if  there  are  multifractal  processes  that  do 
not  possess  the  hallmark  of  SOC  behavior.  Threshold 
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cutoffs  and  finite-size  effects  as  well  as  power-law 
scaling  are  related  to  the  effects  of  different  invariant 
sizes.  Since  ROMA,  by  construction,  can  differentiate 
monofractal  and  multifractal  behavior  in  terms  of  the  in¬ 
variant  scale  sizes,  it  is  a  method  tailor-designed  for  the 
determination  of  this  intriguing  question  related  to  com¬ 
plexity  processes.  The  authors  and  their  collaborators 
have  been  actively  pursuing  the  answer  to  this  impor¬ 
tant  intellectual  question  using  ROMA  along  with  the 
method  of  the  dynamical  renormalization-group. 

5.  Space  plasmas  are  inhomogeneous  at  all  scales.  Ex¬ 
amples  are  found  in  the  auroral  zone  (Uritsky  et  al., 
2008,  2009b),  the  magnetospheric  cusp  as  indicated 
above  and  elsewhere.  As  demonstrated  by  the  calcula¬ 
tions  described  for  the  auroral  electric  field  fluctuations, 
the  rank-ordering  idea  may  be  generalized  to  the  phase 
space  of  multiple  parameters  representing  the  physical 
characteristics  of  the  real  space  environment  such  as  in¬ 
homogeneity.  The  authors  are  currently  applying  such 
an  idea  to  the  cusp  data  observed  by  Cluster. 
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